Exact ground states for a two parameter family of spin 1/2 xyz Heisenberg chains 
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The Heisenberg spin chain (with nearest neighbor interaction) in an external magnetic field, is 
defined by 3 coupling constants (after we re-scale the energy by a multiplicative constant). We show 
that on a particular two dimensional hypersurface, the ground state and all the correlation functions 
can be determined exactly and in compact form. This ground state has a very interesting property: 
all the pairs of spins are equally entangled with each other. In this last respect the results may be of 
interest for engineering long-range entanglement in experimentally realizable finite arrays of qubits. 
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A basic problem in condensed matter physics is to find 
the ground state of a given Hamiltonian embodying the 
interactions of a many body system. A prototype of such 
a system is the Heisenberg spin 1/2 chain, described by 
the following Hamiltonian 



iV 



(1) 



where, cr a 's are the Pauli operators, J Q 's are the coupling 
strengths and B is the external magnetic field. There is 
a long history of attempts for finding exact solutions of 
this model at some special points or lines in this param- 
eter space [1]. The matrix product formalism, originally 
introduced and developed in [2,3], has been recently re- 
vived [4,5] mainly due to the work in quantum informa- 
tion community, where the emphasis is on the properties 
of many body states, like their entanglement or quantum 
correlations [6-8] . In this formalism, one first constructs 
a many body state and then finds the Hamiltonian for 
which this state is an exact ground state. As is always 
the case, when we reverse a difficult problem (in this case 
finding the ground state of an interacting spin system), 
the difficulty shows up in some other form in some other 
place: except for very rare cases, [2] the Hamiltonians 
which arc found are not usually simple and of wide in- 
terest to condensed matter physicists [9]. In this letter, 
we show for the first time that the Heisenberg spin 1/2 
chain can be solved exactly and in compact form on a 
two dimensional surfaces defined by 
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in which (e,?y) = ±(1,±1) are two discrete parameters, 
and g and (J > 0) are two continuous parameters. Un- 
like the xxx or xxz Heisenberg anti-ferromagnetic chains 
whose solutions are implicitly given via the solution of 
Bethe ansatz equations, the ground states of these mod- 
els can be determined quite explicitly and expressed in 
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FIG. 1. In the MPS ground state of (1), all the pairs of 
spins are equally entangled with each other. The curves 
show scaled concurrence for rings of size(from top to bottom) 
6,7,8,9,10,20,30,40,and 50. 



terms of simple functions. Yet as we will see these ground 
states arc quite rich in their properties. We will calculate 
the spin correlation functions exactly and show that sin- 
gularities in the thermodynamic limit develop at g = 0, 
a property which has been called MPS-Quantum Phase 
Transition in [4], to distinguish them from known exam- 
ples of QPT's [10]. 

We will also show that, these ground state have the 
very interesting property that all the pairs of spins have 
equal entanglement with each other. This is a very desir- 
able situation for quantum information processing, both 
theoretically and experimentally, i.e. an array of qubits 
in which there are long range entanglement, figure (1). 
Let us briefly review the MPS formalism. On a ring of 
N sites of d— level particles, a state is called a matrix 
product state if there exist matrices Ai, i = 0, • • • , d — 1 
(of dimension D) such that 
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where Z is a normalization constant given by Z = tr(E ) 



in which E := J2i=o( A * ® A i)- The statc ( 3 ) is re- 
flection symmetric if there exists a matrix II such that 
A[ = IL4jII _1 (where T means transpose) and time- 
reversal invariant if there exists a matrix V such that 
A* = V AiV~ x . All the correlation functions can be cal- 
culated exactly. For example, for a local observable O, 
one finds 
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where E := E^=o(*l b')4* ® A r In thc thermody- 
namic limit (N — > oo), only the eigenvector(s) corre- 
sponding to the eigenvalue X ma x of E with the largest ab- 
solute value matters and any level-crossing in this eigen- 
value leads to a discontinuity in correlation functions. 
Given a matrix product state, the reduced density matrix 
of k adjacent sites is given by 



tr(( A t---A* 
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This density matrix has at least d k —D 2 zero eigenvalues. 
To see this, suppose that we can find complex numbers 



such that 
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This is a system of D 2 equations for d k unknowns which 
has at least d k — D 2 independent solutions. Any such 
solution gives a null eigenvector of p. Thus for the den- 
sity matrix of k adjacent sites to have a null space, it is 
sufficient (but not necessary) that d k > D 2 . Let the 
null space of the reduced density matrix be spanned by 
the orthogonal vectors \e a ),a = l,---,s, then we can 
construct the local Hamiltonian acting on k consecutive 
sites as 



a=l 



where J a are positive constants. The total Hamiltonian 
on the chain will then be given by thc positive opera- 
tor H = '}2d =1 hijj r k, where /iu+fc is the embedding 
of h into sites / to I + k of the chain. The state \ip) will 
then be a ground state of H. Equation (d k > D 2 ) puts a 
stringent requirement on thc dimensions of the matrices 
used in construction of a matrix product state. When 
dealing with spin 1/2 with nearest-neighbor interactions, 
for which d = 2 and k = 1, it appears that the only ad- 
missible dimension for the matrices Ao and A\ is D = 1, 
leading to a product state. However, it is crucial to note 
that the condition d 2 > D 2 is only a sufficient and not 
a necessary condition for the density matrix p to have a 
null space. 



To proceed with our construction, we require that the 
state satisfy some natural symmetries, i.e. spin-flip sym- 
metry which, in the language of matrix product formal- 
ism, means that there is a matrix X such that 



XAnX- 



eA lt XAiX- 1 =eA , 



where e 2 = 1. Working in the basis where X = er 2 , we 
find the general form of the matrices Aq and A\\ 
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Although these two matrices are not symmetric, the state 
constructed from them is symmetric under parity, since 

with thc property 



there is a matrix n 



iL4*rr 1 = Ar 







rL^ir 1 = Al 



We now consider the matrix equation (5) which in the 
present case is 



cqoA) + cqiAqAx + aoAxAo + cnA\ = 0. 



(6) 



This is a set of linear equations for the four coefficients 
which can be written as a matrix equation MC = 0, 



leading to a non-zero solution when 

det(M) = 16b 2 c 2 (a - d) 2 (a + df 



0. 



Thus we will find non-trivial models, for a = d or a = —d. 
The models with b = or c = are not symmetric un- 
der parity, since in these cases the matrix LT will not be 
invertible. We can always re-scale the matrices by a con- 
stant factor without affecting the matrix product state, 
so we set a = 1 and use a subsequent gauge transfor- 

c 
1 

Therefore we are left with the following four classes of 
models defined by the matrices 
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where g is a continuous parameter and (e, rj) = ±(1, ±1). 
The four types of models arc distinguished by the val- 
ues of the pair (e, 77). The eigenvalues of the matrix 
E = A (g) A + At <g) Ai are 2(?y ± g), 2(1 ± g) . The 
correlation functions can be derived from (4). Consider 
for definitcness, the case 77 = 1. The magnetization per 
site is found from (4) to be 



K) = (<j z ) = 0, 
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where u := jt^- and the correlation functions G a (l, r) 
{<Ta,iGa,r) ar e similarly found to be as follows: 
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FIG. 2. (Color Online)The magnetization in the x direction 
as a function of g, (e = 1). 



G x (l,r) = 
G y (l,r) = 
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These correlation functions satisfy the following rela- 
tions: 

G x + G y + G z = l , (1-G z )(l-G y ) = (a x ) 2 . (9) 

In the thermodynamic limit (N — ► oo), discontinuities 
develop in these correlation functions at g = 0. For ex- 
ample, the magnetization per site will be 
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and displayed in figure (2). Interestingly the spins align 
themselves opposite to the magnetic field. To under- 
stand this, let us set for definiteness e = ry = 1 and 
consider the two-body Hamiltonian i?2, after subtracting 
the ferromagnetic term —Jo ■ a for which any product 
state \<t>)\4>) is an eigenstate. For \g\ ^> 1 the remain- 
ing Hamiltonian tends to H2 ~ g 2 {\&ix&2x + <?ix + &2x) 
, with the ground state \x— )\x— ) and for \g\ 1 it 
tends to H2 ~ {^&ix&2x — <*\x — &2x)> with the ground 
state It is thus the combination of the anti- 

ferromagnetic interaction in the x direction and the mag- 
netic field which align the spins opposite to the magnetic 
field. Even when we tune the coupling J so that there is 
no xx coupling, (i.e. J = (1 +g 2 )/2 ) this ant i- alignment 
happens. To understand this, consider g w 1, where 
H2 ~ — &iz a 2z + B{o\ x + <J2x) with B very small. The 
ground state of this Hamiltonian is \ip) « \4>+) ~ B\ip + ), 
where \ip+) and \4>+) arc Bell states. One then finds that 
{(Tlx + o~2 X ) = —^B, confirming figure (2). 

In order to see how the Hamiltonian is constructed, we 
solve equations (6) which in view of (7) take the form 

(1 + g)(C 00 + C n ) + e(l - <?)(C i + Go) = (10) 
(1 + V )(C 00 - C n ) - e(l - r?)(C i - C 10 ) = 0. (11) 

It is easy to verify that the solution space is determined 
by the following two un-normalized vectors, 



\e 1 ) = (l+ V )\i>-) + (l-r,)\4>-) 

\e 2 ) = (l+gM+)-e(l-g)\cb + ), 

where |V>±) = ^(|0,1) ± |1,0)) and |0±) - ^(|0,0) ± 
|1,1)) are Bell states. Under spin flip the above states 
transform as |ei 2) — * T-|ei.2)- The final local Hamilto- 
nian will be given by h = J|e 1 )(e 1 | + |e2) (e2 1 , where J is 
a non-negative parameter and we have used the freedom 
for rescaling the couplings of the Hamiltonian to set one 
of the parameters equal to 1. In view of the symmetry 
property of the vectors under spin flip, this Hamiltonian 
will be symmetric under spin flip. Expressing the above 
operator in terms of Pauli operators and subtracting con- 
stant terms, we find the total Hamiltonian which is writ- 
ten in (1) and (2). 

Note that the models with e = ±1 (sign of the magnetic 
field) are related by local 7r-rotations of spins around the 

Also, the models with 



axis where (cr 



■r.y 



'x,yl 



77 = ± are related to each other by simultaneous rota- 
tions i?x(f ) <£> Rx(-^-) of spins on adjacent sites, under 
which we have {<J z ^u z .i+i ~ cr y,iVy,i+i)- This is, of 
course, possible only when N is even [13]. 

The explicit form of such a ground state can also be 
determined. For rj — 1, the matrices Aq and A\ commute. 
By a similarity transformation which does not change the 
state (3), both the matrices are made diagonal, 

A °-{ l-^J'^'K 1 + V g) 



and the MPS state (3) will be given by 

'*>'=' " 71 

where 



(12) 



|^±) = (i±V?)|o) + (itV5)|i), 

and Z = 2 N+1 ((1 + g) N + (1 - g) N ). These expressions 
are valid for all values of <?, provided that we replace 
Jg — > i-\f—g when we consider negative values of g. 
Note that (</>+|0_) = 2(1 — g). One can indeed check that 
the separate product states are ground states of H, (i.e. 
the local Hamiltonian acting on two adjacent sites, when 
added by a suitable constant, annihilates \4>±)® 2 )- How- 
ever the advantage of the MPS state we have constructed 
10+}®* + \<f>-) 9N or the other state \(f> + )® N - 1^-}®* 
is that they are invariant under spin-flip transformation 
(X® . Thus even if the couplings of the Hamiltonian are 
not tuned as exactly as in (2), and are perturbed a little 
bit, first order perturbation theory guarantees that one 
of these entangled states and not the product states, will 
be the unique grounds state of H. 



We now come to the entanglement properties of the 
state (12). At g = 1 when = 0, the state be- 

comes a standard GHZ state, ^~(|0 • ■ ■ 0) + 11 • ■ ■ 1}). For 
other values of <?, when \<p+) and \<f>—) are no longer or- 
thogonal, it can be named a generalized GHZ state. Ob- 
viously such a state induces equal entanglement between 
any two spins regardless of their distance. To calculate 
this entanglement we determine the reduced two particle 
density matrix and use Wootters formula [11], with the 
result [12]: 



C 



4|ff| 



l + g)N + (l-g)N\ 
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Thus although the ring is not totally connected, the mu- 
tual entanglement of all pairs are equal and independent 
of their distances. Looking at the N 3> 1 limit, one can 
obtain the relation 



v iV' ' cosh g 



I si 



(13) 



One can interpret the left hand side as the total mutual 
entanglement of a spin with all the other spins, and the 
above equation as a universal scaling relation for this to- 
tal entanglement. 

For the case r) = — 1, we use the transformation 
U : = i? K (f ) <£> Rx{^-) on adjacent sites which trans- 
forms H(rj = 1) to H{r) = -1) and act on (12) by U® 1 ? 
to obtain the ground state as 



where 



1 



{\x+)\x-)f^ + {\x-)\x + )Y 



(14) 



|x±> = (l + V?)|?/,±) ± i(l-^g)\y,T), 

in which |y,±) denote the eigenstates of <r y . An alter- 
native way for deriving this ground state and indeed the 
reason for its simple structure, is to note that for r\ = — 1 , 
although the matrices Aq and A\ do not commute, the 
pairs of matrices corresponding to Bell states, defined by 



V2 



(A A m + (-l) n A 1 A 1+m ), m,n = 0,1 



commute with each other. The reader can verify that the 
states |x)±|x=f) are indeed linear combinations of the Bell 
states (poo = 4>+ , (poi = ip- and <pio = V J +i making the 
state (14) a linear superposition of strings of various Bell 
states on adjacent sites. 

In summary we have introduced a two-parameter family 
of spin 1/2 xyz Heisenberg chains, with nearest neigh- 
bor interactions in an external magnetic field, for which 
ground states and all correlation functions can be cal- 
culated exactly. These states have two very interesting 



properties: first, they undergo a discontinuous (quan- 
tum phase) transition as one of the parameters passes 
a critical point, which stimulates further exploration of 
MPS-quantum phase transition in a set of important ex- 
actly solvable models. Second, they have the property 
that all the pairs of spins are equally entangled with 
each other. This makes them good candidates for en- 
gineering long-range entanglement in experimentally re- 
alizable arrays of qubits or spin systems. This study can 
be extended in several directions, including generaliza- 
tion to open chains, finding the excitations, perturbing 
them to explore more models, and actual engineering of 
such chains for small number of qubits for information 
processing tasks. We thank David Gross, of Imperial col- 
lege, London for a very stimulating email correspondence 
which led to substantial improvement of this paper and 
also A. Langari and M. R. Rahimitabar, for their valu- 
able comments. Corresponding author, V. Karimipour, 
vahid@sharif.edu. 
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